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Matrix A: Standard Mass Matrix, Easy 


Notes: 
‚ш B: Usually Requires Union Mesh 


. FOR Non-Conformal Meshes: NUMERICAL RESULTS 
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R.O.S.E. FOR DDMS WITH CURVED INTERFACES 





Original TCDA unit cell 
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Original/Planar Array Features: 


> Tightly coupled dipole array with 
FSS superstrate is shown on above 


» Operational Bandwidth : 
0.5-3.25GHz 


> Max. scan-angle +70° in all planes 
normal to aperture 














Antenna design taken from Prof. Volakis’s 
group 








Finite array on planar ground plane 





Unit Cell Conversion: 
Mapping from Planar to Cylindrical Surfaces 


PGS sees 


J › LE 
х-рсовф,у-рвіпф (х „У 2) 










R.O.S.E. FOR DDMS WITH CURVED S 
INTERFACES | 








Note: Difficult to construct Union mesh for non-conformal 
meshes on curved surfaces 











» Run with FEM-DDM R18, IEDG EXT R11 
> Running frequency f = 2 GHz, curvature radius r = 1500mm 


1 


IE region mesh 


diation pattern [dB] 
ж ROSE ТЕТІ 
rmal 


Model 1 — Conformal with FEM Model 2 — Coarse mesh (A/15) 
outer surface mesh compared with FEM outer surface 
mesh (ranges from A/25 to A/15 ) 


s Modi — IMod2 
Solver MS-DDM r31 ROSE 
Iter. 8 (conv. 0.0098) 7 (conv. 0.0063) 
Wall time 14h27m46s 15h45m13s 
RAM 50GB 19GB 








CEM TECHNIQUE: BOUNDARY ELEMENT METHOD (MOM) 


Representation Formula 
E*(x) 5 jkn [ gx. y)J(y) dy - 7v | g(x, y)V;J(y)dy, x ex 
до до 
Taking the limit x — x e 0Q, we have 
z E° (x) Š A (J) 
A = JT] s ^ 
= јкпл, | о(њуј(у)ду + V, | eoe y) V Ji) dy 
до до 


with IV=nxvxn 


EM Scattering: Assuming 20 is РЕС 






Define a residual magnetic current М(х) 
M(x) = nx E™(x)- n x A (J), 
Note that M(x) е H(div,;T) 


ue H(div,.,T) = № (x) =e” (x) — oin, (x) + VA, E) 
should be tested by ALL possible currents 
=> Find ue H(div,.,T) such that (V, RE ). - 0,Уу е н(фуГ) 





|| |выв :тЕ"(х)= ката) Буто) 
Find ue H(div T) such that 


h (v.e ). = ж Қа, )- 7 (ç, vo, 3 
Му е H(div T) 











U Finite Dimensional Discretization 


~ _ zinc я : : 
Zort =E", Zr: An NxN complex, symmetric, and dense matrix 


Jin-Fa Lee Lecture Notes 













IE-DG 
SIMULATION AT X- 
BAND, 

10 GHZ, USING 
3GHZ MESH 


1-to-9 automatic 
h-refinement 


IEDG FORMULATION FOR NON-PENETRABLE TARGETS 





16.4.5 Galerkin Weak Formulation 


We now have a linear combination of the weighted residuals as: 


Cy GRO) + са (v, RO) + 8 (é : v, RG) + сз (4 Е v, RO) = 0 
i I (16.42) 





which includes both the EFIE and MFIE surface penalty and the boundary 
penalty terms. In general, the choices of c4, со, 8 and cs can be determined by 
the considerations of accuracy and/or convenience. Similar to the conventional 


СЕТЕ formulation for PEC non-penetrable objects, herein we set cq = c2 = i. 


Our choice of c3 is mainly the consequence of convenience. We have chosen 
сз = -4 such that the problematic contour integrals, appearing in the fourth 
inner product term of Equation 16.27, can be canceled out. Moreover, the 
interior penalty stabilization function 8 is taken as 8 = ah~', where h is the 
average mesh size of the discretization and the stabilization parameter a is a 


positive number independent of h. Throughout the numerical examples reported 
in this paper, we have employed a = т and subsequently, we have the interior 


penalty stabilization function set to be 5 — ih |. 





Notes: a. Still integration by part and need contour integrals; b. The contour penalty terms are un- 
orthodox and difficult to explain (although absolutely correct) 








Note: Unbounded Domain + Silver-Muller Radiation Condition => UNIQUE 


E,H 
n 


~ Boundary Value Problem (a): 
За s n s S 
V x H!) = J@eE + J i У х НО = уске 


) 
\7 / lim ЦУхЕ jx xt?]- 0 


(ко но) Boundary Value Problem (b): 
2 
VxEP--jouH ` УХЕ = - jguH 
I , m О, апа 
У х НО = јосЕ V x HS” = jock!” 


in Q 


2 











lim Ix (Vx? xx Et" 20 


je 











ге 
Е 
| 
= 
г 7 
e 
=> 


N | nx(H-H”)=J 
5 


a оп Г 
~ м” (Е?-Е”)ха-м, 








SURFACE EQUIVALENT THEOREMS 


Q: How to set for J, and M, such that (ЕФ, но ) = (0,0) апа (E^. nt (EY), HS”)? 


js 
J, -nx(H? - H)- nx(Ht?- 0)= nx H 
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REPRESENTATION FORMULAE IN FREE SPACE 
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FINITE DISCRETIZATION: BEM IMPLEMENTATION 











FINITE DISCRETIZATION: BEM IMPLEMENTATION 
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R.O.S.E. ALGORITHM FOR COMPUTING SELF-TERM 
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R.O.S.E. ALGORITHM: CONT 
From Equivalent Principle => Three Linearly Independent Plane Waves: [E^ „но | ДЕФ HO} {E®, HO} 
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R.O.S.E.: Numerical Results 


В-о (RWG Basis Functions) 
Interpolation 


Notes: 
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Figure 4: RCS results of the sphere: (a) RCS comparison of different solvers; (b) relative L2 


error convergence 











R.O.S.E.: Numerical Results 






J Magnitude Table 2: Iteration time versus electrical size for the Cube 
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R.O.S.E.: Numerical Results 
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R.O.S.E.: Choices of Basis Functions 


Notes: 
1. EFIE with traditional integration by part 
2. 8B =a for both EFIE-inter and EFIE-proj 
3. B = BC Basis Functions for EFIE-mix-proj 











Notes: near field RMS errors for two different frequencies for a PEC sphere 





RMS error 
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R.O.S.E.: Error Investigations 










Given an current distribution J, on I’, the reaction of J, with a test current 
J, can be defined as 
lea = / J; e E'^*(J ,)dr 
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Figure 5: The relative L2 error convergence of the self terms and reaction terms of ROSE 


2 š 
L^ relative error of current 








R.O.S.E. for EFIE: no integration by part 
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R.O.S.E. with EFIE: no integration by part, self and touching terms 





Unit sphere at 150 MHz 





(a) 


Fig. 4. Singularity: (a) self-term singularity, basis RWGs share at least one 
patch with test RWG, NV; = 5; (b) near singularity, basis RW Gs share at least 
one edge with test RWG, NV; = 11, 12, 13. 
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R.O.S.E. with CFIE: no integration by part, self and touching terms 
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